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We prove a topological criterion for the existence of zero-energy Majorana bound-state on a 
disclination, a rotation symmetry breaking point defect, in 4-fold symmetric topological crystalline 
superconductors (TCS). We first establish a complete topological classification of TCS using the 
Chern invariant and three integral rotation invariants. By analytically and numerically studying 
disclinations, we algebraically deduce a Z2-index that identifies the parity of the number of Majorana 
zero-modes at a disclination. Surprisingly, we also find weakly-protected Majorana fermions bound 
at the corners of superconductors with trivial Chern and weak invariants. 
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Symmetry protected topological insulators and su- 
perconductors have theoretically, and experimentally, 
risen to prominence in the last half-decade pQ. Re- 
cent developments in this field have moved on from the 
study of discrete symmetries such as time-reversal and 
charge-conjugation|2 -;4j, to translational and point group 
symmetries [SHHj. While spatial symmetries are not pre- 
served as generically as, say, time-reversal, they can still 
support robust topological states in clean crystalline sys- 
tems. It is understood that the so-called strong topo- 
logical invariants, which are protected no-matter what 
spatial symmetries are broken, determine the appearance 
of disorder-insensitive gapless boundary states. Interest- 
ingly, it was found that weak invariants, which require 
an additional translation symmetry, support boundary 
states [TSJ [TB], and more surprisingly, robust bound-states 
on crystal dislocations [PTrl2"l"] ■ A natural question is then 
to ask whether topological defects of the point-group ro- 
tational symmetry, i.e. disclinations, can also bind low- 
energy states in topological phases protected by point- 
group symmetry. A related problem has been studied in 
a different context in graphene |22H24j . In this Letter we 
address this question for 2d topological superconductors 
with point-group symmetry. To determine the existence 
of zero-energy Majorana bound states (MBS)[25j ES] on 
disclinations we derive an index theorem connecting the 
parity of the number of zero modes to the eigenvalues of 
rotation operators at rotationally invariant points in the 
Brillouin zone (BZ). Our derivation combines a unique 
algebraic approach with exact diagonalization numerics 
and can be generalized to prove generic index theorems. 

Just as dislocations are local topological defects in the 
translational order of crystals, disclinations are topolog- 
ical defects in the discrete rotational order. In 2d, discli- 
nations are point defects that can be constructed by a 
Volterra process [37] of removing or inserting material in 
certain angular sections with angles compatible with the 
crystalline symmetry (see Fig. [T]) . Dislocations are char- 
acterized by their Burgers' vector i.e. the translation 
element acquired when a particle encircles the defect; 
disclinations are described by an element in the space 
group encoding the amount of rotation Q and transla- 




FIG. 1. (a), (b)Lattice configurations with three disclinations. 
Each of the three faces is a L x L square lattice. Periodic 
boundary conditions are taken on the six edges of the cubes 
as indicated by matching colors/line styles, (c)-(f) Flattened, 
zoomed-in versions of the —90° (c,d) and +180° (e,f) discli- 
nations at O, K and K' . (c) type-1 disclination centered at 
a triangular plaquette with an odd number of translations 
around its boundary, (d) type-0 disclination centered at a 
trivalent vertex with an even number. 

tion T one picks up traveling around the point defect. 
The translation piece T is not unique and depends on 
the enclosing path, however, for a C4 symmetric lattice 
the evenness (type-0) or oddness (type-1) of the number 
of translations is unique. The C4 symmetry thus yields 
a Z2-characteristic that distinguishes classical disclina- 
tions with the same Frank angle Q. Examples are shown 
for —90° disclinations in Fig. [TJ:,d. The difference be- 
tween the two primitive disclinations can be observed at 
the defect cores where C4-symmetry is violated at either 
a triangular plaquette or a trivalent vertex. The Z2 char- 
acteristic arises from the fact that there are two inequiv- 
alent 4-fold rotation centers, vertex or plaquette, which 
can be extended to general C„ [2"5], 

We will only consider fully-gapped, translationally in- 
variant superconductors in the mean-field limit which 
are described by Bogoliubov-de Gennes (BdG) Hamil- 
tonians H(k), in Bloch form, with a particle hole con- 



straint Si/(k)S^ = — H{— k), where S is a local, anti- 
unitary operator. A point group element r is repre- 
sented by a unitary operator f that commutes with the 
full Hamiltonian, and satisfies rH(k)f^ = H(r'k) for the 
Bloch Hamiltonian, and SfE^ = f for the particle-hole 
operator [22]. For the duration of our work we will fo- 
cus on the Abelian point-group C4, generated by ir/2- 
rotations, which is a symmetry commonly shared by all 
layered perovskites. We use the half-integer spin rotation 
such that f 4 = (f 2 ) 2 = — 1. 

^-symmetric superconductors in two dimensions is 
classified by (i) the Chern invariant 
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(k)|<9 fe Ju„(k)) being the Berry 
connection of the negative-energy bands [30 32 , (ii) the 
eigenvalues of the rotation operator f for all the neg- 
ative energy states at the 4-fold symmetric momenta 
IT = r, M(= (71", 7r)) (iii) the spectrum of the C2 rota- 
tion f 2 at one of the equivalent 2-fold symmetric mo- 
menta X = (w, 0) or X' = (0,7r). We label the bands 
at the symmetry points by their rotation eigenvalues 
following Ref. 133] At II = T,M, a band with f = 
e -W4 ( e «/4 ) e 3W4 ) e -3i7r/4 ig labelled by n 5 ,n 6 ,n 7 ,n 8 



respectively; while at X, a band with r 



is la- 



belled by X 3 , X4 respectively. Let #Hi, #Xi be the num- 
ber of appearances of IL;, Xi within the negative-energy 
states. Only differences of the eigenvalues between sym- 
metry points carry topological information. Our conven- 
tion uses the eigenvalues relative to the values at the 
F-point, and we define: 



n 3 = #A 3 
714 = #A 4 

m = #m, 



#F 6 - 
#F 5 - 
#F 2 , 



#F 8 
#F 7 
for i 



5,6,7,8. 



(2) 
(3) 
(4) 



These are easy to understand as r 6 ,r 8 (r 5 ,r 7 ) both 
square to +z(— i) and thus 773,774 determine the differ- 
ence in the C'2 eigenvalues at X and T while 775,6,7,8 de- 
termine the C4 eigenvalue differences between M and 
r. Including ch this gives seven topological numbers, 
but there are constraints which imply the topology is 
determined by fewer quantities. These integers obey 
71 3 + 774 = n 5 + uq + nj + n$ = (from the constant num- 
ber of negative-energy bands throughout the BZ) and 
715 + 7i6 = 71 7 + 7i 8 = (from the particle-hole constraint 
ensuring the rotation spectrum at symmetry points of 
the unoccupied bands is the complex conjugate of that 
of the occupied ones). Following the work of Ref. [T4lone 
can show that 



ch + riQ + 2n4 + 37i7 = mod 4. 



(5) 



Hence C4-symmetry topological crystalline superconduc- 
tors are completely classified by four integral invariants 



FIG. 2. C4-symmetric tight binding models with four Majo- 
rana fermions (black dots) at each site, (a) Hb (b) H c 
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TABLE I. Chern and rotation invariants for the chiral p x +ip y 
superconductor |7f . The Hamiltonians are superscript labeled 
by their Chern and weak invariants (ch; u) . 
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TABLE II. Chern and rotation invariants of models in Fig. [2] 

Xi = (ch, 714, n 6 , 717) that satisfy ([5]). Moreover the weak 
Z2-topological invariant is determined by the inversion 
eigenvalues at M and X: 

G„ = f(Gi + G 2 ), v = (774 + 7t 6 - 717) mod 2 (6) 

where Gi, G2 are the reciprocal lattice vectors. [51] 

The essential ingredient of our index-theorem proof is 
a collection of 2d C4 symmetric superconductor models 
that "generate" all the topological classes characterized 
by the Xi- Since Xi is a f° ur component vector we need 
four Hamiltonians which have linearly independent Xi- 
Combining the Hamiltonians via direct sum combines the 
vectors with a vector sum, so any topological class x% can 
be produced. We choose two generators to be spinless, 
chiral p x + ip y superconductors on a square lattice: 

H a = A (sin k x r x + sink y T y ) 

+ [iii(cos k x + cos k y ) + 2u 2 cosk x cos k y ]r z (7) 



where r a acts on Nambu-space, A is the p x + ip y pairing 
and the first/second neighbor hoppings Ui,U2 are kinetic 
energy terms that gap out the nodes of the pairing term 
at the points T, X, X' , M. The particle- hole and rotation 
operators are given by S = t x K and f a = (I2 +7Tz)/v / 2 
where K is complex conjugation. The invariants \i de- 
pend on Tii and 7i2 and are summarized in Table [I] Flip- 
ping the signs of both U\ and 7i 2 inverts \% ~ > ~ Xi- 

The other independent generators are 2d general- 
izations of Kitaev's p-wave wire [351. Fig. [2] depicts 



two tight-binding limits of C4 symmetric Majorana 
fermion models with four fermions per site and arrows 
which represent the Majorana ordering convention. The 
two Hamiltonians are H b = ^E x (7x7 x + ei + 7x7 x +e 2 ) 

and H e = itJ2 x (7x7x+e 1+e2 +7x7x-e 1+e2 )> where 

7x's are Majorana operators with 7^ = 7^ and 
{7x)7y} = 2<J v <5 xy . The C4 rotation operator 
he = rix ex P(-f7x7rx)exp("f7x7rx)exp(-f7x7rx) 

gives hcilL 7x>7 X )7x)C = (7rx: 7r X , 7rx, -7, 4 X )> where 
r is the C4 rotation in real space. If we transform to 
complex fermions c = (7 1 + ?7 3 )/2 and d = (j 2 + ij 4 )/2 
then the corresponding BdG Hamiltonians take the block 
diagonal forms 

iJb(k) = t(cos k x r z + sin k x T y ) © t(cos k y r z + sin k y T y ) 
H c (k) = t(cos(k x + k y )r z + sin(k x + k y )r y ) 

ffi t(cos(k x - k y )T z + sin(fc y - k x )r y ) (8) 

in the basis £k = (c_k, c^, d-k, d\ l ) T and where r^'s 
act on the Nambu degree of freedom. The particle- 
hole and rotation operators are 5{, c = (I2 ® t x )K and 
He = 0"+ © I2 — io~- © Tz where <j± = l/2(a x ± ic,,) 
acts on the (c, d) space. The lattice structure of H c is 
two copies of Hb displaced by half a lattice spacing, how- 
ever, as a ^-symmetric electronic structure, H c (k) is 
not two copies of Hb(k) because of the different rota- 
tion centers. This is evident in the respective \i shown 
in Table [TTJ Together with the two tight binding models 
H^ 1 ' ^ and Ha' 1 ', Hb and H c generate all possible com- 
binations of Xi — { c h,n^,nQ,n-j) allowed by |5]). Thus, 
every ^-symmetric superconductor must have the iden- 
tical topological properties to direct sums of the two H a 's, 
Hb, and H c . In particular the zero- modes at disclinations 
of any ^-symmetric superconductor can be determined 
from these four models. 

After constructing the Hamiltonian generators, we 
must determine the properties of disclinations for each 
generator. A disclination configuration is specified by the 
pair (f2, T) where is the Frank angle and T is the trans- 
lation. The MBS of the chiral superconductors H a (with 
pairing and hopping parameters 2^2/ A = ±ui/A = 1) 
are studied numerically using periodic lattice models 
with three disclinations. This is achieved by taking three 
adjacent faces of a cube and gluing the parallel sides 
(see Fig. [lji,b). Two —90° disclinations are located at 
the points O and K of positive curvature and a +180° 
disclination is located at the point K' of negative cur- 
vature. We use two lattice configurations that differ in 
the choice of a type-0 or type-1 disclination at O (see 
Figjl]:,d). The superconducting phase is chosen so that 
it smoothly winds around the defects, but there are two 
inequivalent ways of specifying the defect lattice due to 
the double covering of the rotation group. The smooth 
winding around the disclination involves either the ro- 
tation f(s) = e ls ° r z/2 or ^S'( s ) — e Js(27r+n)r 2 /2^ p arame _ 

terized by s € [0, 1]. We choose f(s) and thus the phase 
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FIG. 3. (a,b) Exponentially localized probability amplitudes 
of the Majorana zero mode at disclinations O and K, plotted 
on a torus geometry where parallel sides on the hexagonal 
domain are identified. Compare with 3d lattice in Fig. [T] 
(c) Tight binding model Hb with (left) type-1 and (right) 
type-0 disclinations (d) H c with (left) type-1 (right) type-0 
disclinations. Thick red dots in disclination cores are unpaired 
Majorana boundstates. 

smoothly winds by 7r/2 around O, K and — 7r around K' . 
Note, if we choose f'(s) then the corresponding 4- fold op- 
erator would be — f instead and the invariants are flipped, 
rig <-> —77.7; the index is then calculated the same way af- 
ter this change. We find that only H a X supports an 
odd number of MBS and even then only for type-1 discli- 
nations. The results are summarized in Table ITTT1 and we 
show the zero mode wave functions at the disclinations 
at O and K in Fig. |3^,b. 

For Hb, H c we use the fact that the parity of the num- 
ber of MBS at the defect is insensitive to all perturbations 
that do not violate the energy gap or rotation symmetry 
away from the point defect since there is no low-energy 
channel for a single MBS to escape or enter the discli- 
nation core. This implies that, just as for the boundary 
modes of the topological p-wave wire [35 , we can deter- 
mine the parity of the zero-mode bound states pictorially 
in the tightly-bound limit. The MBSs of Hb and H c at 
type-0 and type-1 —90° disclinations are shown in Fig. 
[3];,d, where the MBSs are simply un-bonded Majorana 
fermions represented by thick red dots. We find that Hb 
has a zero-mode for type-0, and H c has zero modes for 
both types. This is summarized in Tabic |Hl} 

We are now in a position to determine the index the- 
orem since any C4 symmetric BdG hamiltonian can be 
smoothly deformed into a unique composition 

[H] ~ mi [H^] © m 2 [H^} © m 3 [H b ] © m A [H c ] (9) 
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TABLE III. Parity of the number of zero modes at a 
disclination for the chiral superconductors H^~'°\ Hi 1 ' 1 ' in 
with smooth rotation f(s) 
models Hb,H c in Fig. [2] 



90° 


jsflTt/H tight bindmg 



up to topologically trivial bands far from Fermi-level. 
The nii are integers, and m, < means a direct sum 
with the negative and positive energy states switched. 
To finish the derivation of the topological index we need 
to carry out some simple algebraic manipulations. First, 
we see that since H c only has 714 = 2 non-zero and has 
zero-modes for both disclination types we determine that 
the index gets a contribution of 1/2(714) mod 2. Next we 
can take 2i^ 1;0) © 2H£' A) © H c which has Xi = (4,0,0,0) 
and bound-states for both types. This implies the in- 
dex recieves a contribution of l/4(c/i) mod 2. Using these 
two pieces we can go back to Hi l '°\ which does not have 
any zero-modes, and determine the equation [1/2(714) + 
l/4(ch)+k(rie)] mod 2 = which upon substitution gives 
k = 1/4. To determine the contribution of 777 we consider 
H { a m © © H b with X i = (2,0,-1, 1). This model 

has bound states on both types of disclinations so we 
use [1/2(714) + (l/4)(cft) + (l/4)n 6 + jn 7 ] mod 2 = 1 
to find j = 3/4. So far we were careful to choose 
all the Hamiltonian combinations above to have a van- 
ishing weak invariant as it also contributes to the in- 
dex. We can see this by taking H^ 1 ' 1 ^ which yields 
1/4 [c/i + rc.6 + 2ri4 + 3-127] mod 2 = 0, yet has a bound 
state on type-1 disclinations. This bound state, however 
arises from a different mechanism, and comes from the in- 
terplay of the non-zero weak invariant and the oddness of 
the translation T around a type-1 disclination. Thus, we 
have determined the existence conditions for an odd num- 
ber of MBSs at a —90° disclination. Combining disclina- 
tions gives (Q 1 ,T 1 ) + (Q 2 ,T 2 ) = (ni+n 2 ,Ti+r(fii)T a ) 
which implies that for generic C4 disclinations with Frank 
angle f2 the topological index is: 



= 



1 T • G v + — (ch + n 6 + 2n A + 3n 7 ) 



mod 2 



2tt 2tt 

(10) 

where G„ is the weak Z 2 invariant. As mentioned, the 
first term resembles the topological index for MBSs at a 
dislocation, where the Burgers' vector B is replaced by 
T [rTWT§] . It vanishes for type-0 disclinations and equals 



714+77-6+717 (mod 2) for type-l's. The second term of (10) 
is an integer because of the constraint ^ on the Chern 
invariant and can distinguish Chern numbers which are 
even or odd multiples of four. 

For a more physical understanding we can refer to the 
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FIG. 4. Untwisting a —90° disclination of the four layer p x + 
ip y continuum model Hi into a single layer p x + ip y model 
with a quantum vortex. 



outer boundaries of the regions surrounding the discli- 
nations in Fig. [3j;,d. For Fig. [3j; we see the boundary 
links carry one unbound Majorana and the corners con- 
tain two whereas in Fig. [3]i the links carry two and the 
corners carry three. Since there is only one disclination, 
the parity of MBSs on the boundary will match the par- 
ity in the disclination core. We can clearly see that the 
two terms that comprise represent the edge and cor- 
ner contributions to the index respectively. T • G u /2it 
counts the number of MBSs (mod 2) on an edge with 
length T, while (ch + ?76 + 2774 + 3777) /4 counts the num- 
ber of MBSs at a 90° corner. The index ( 10 ) therefore 



not only gives information about the disclination core, 
but also the defect-free system boundary. In particu- 
lar, even a system such as H c , with vanishing Chern and 
weak Z 2 invariants (which thus does not carry topologi- 
cally protected edge modes) binds Majorana fermions at 
corners since = 1. This implies the existence of MBSs 
in the form of corner states even in a defect and vortex 
free system. 

can be illustrated in a continuum system on a disk 
geometry with a disclination at the origin. Take four 
copies of a spinless, continuum p x + ip y superconductor 
with each copy rotated by 7r/2 relative to the previous: 



V/2 



H4 — ho 
fto(k) = 

discrete 4-fold rotation symmetry 



A|fc^T r + \A\kyTy + (777 



for h$ = e^/ 2 /7 e-^/ 2 , 
ek 2 )T z . H 4 has the 



( l a \ 

a 2 

i 2 

V-12 / 



(ii) 



A corresponding —90° disclination is represented by the 
helix in Fig. [4j where the top and bottom layers are 
glued along the branch cut (red line) with anti-periodic 
boundary conditions since f| = — 1. The disclination he- 
lix can be untwisted to form a single copy of a p x + ip y 
model with a 7r-flux vortex replacing the disclination at 
the origin that binds a MBS|25j. This is consistent with 
the index theorem ( 10 1 since ch — 4, 774 = T = for con- 
tinuum models, and also 77 6 = 777 = since momentum 
space can be compactified into a sphere S 2 = R 2 U {00} 
where the rotation spectra at the fixed points k = and 
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oo are identical. 

We have shown that topological crystalline supercon- 
ductors in two dimensions with C4 rotation symmetry 
are classified by four integers \i = (c/i, J24, ri6, 717). The 
appearance of zero-energy MBSs at disclinations is deter- 
mined by a Z2 topological index 0. Although the index 
theorem relies on C4 symmetry, a MBS is robust against 
any rotation breaking perturbation that does not close 
the bulk energy gap around the defect. In fact, the discli- 
nation itself breaks rotation symmetry. The MBS can be 
interpreted as a consequence of the topological nature 
of the BdG hamiltonian with a defect -ff(k,s)[TH] which 
does not assume any rotation symmetry (here s adiabat- 
ically parametrizes a loop around the point defect). The 
MBSs at the corners of a sample, however, are sensitive 
to C4 breaking perturbations, and can escape through 
accidental low-energy edge channels. If the edge recon- 
struction is weak enough, MBSs at corners could provide 
the possibility for indirect observation of the elusive ex- 
citation through transport experiments. 
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